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ABSTRACT 

A tidal disruption event (TDE) takes place when a star passes near enough to a massive black hole to 
be disrupted. About half the star’s matter is given elliptical trajectories with large apocenter distances, 
the other half is unbound. To “circularize”, i.e., to form an accretion flow, the bound matter must 
lose a significant amount of energy, with the actual amount depending on the characteristic scale of 
the flow measured in u nits of the black hole’ s gravitational radius (~ 10®^(i?/I000i?g)“^ erg). Recent 


numerical simulations (Shiokawa et al. have revealed that the circularization scale is close to the 


scale of the most-bound initial orbits, ~ \ 5 ^g ~ cm from the black hole, and 


^BH,6.5 
^44 


the corresponding circularization energy dissipation rate is ^ 5 erg/s. We suggest that the 

energy liberated during circularization, rather then energy liberated by accretion onto the black hole, 
powers the observed optical TDE candidates. The observ ed rise times, lumin osities, temperatures. 


emission radii, and line widths seen in these TDEs (e.g. Arcavi et al. 20141 are all more readily 


explained in terms of heating associated with circularizatioh than in terins of accretion. 


1. INTRODUCTION 

Occasionally, a star passes close enough to a super- 
massive black hole (SMBH) to be tidally disrupted. The 
frequency of such occasions is estimated to be ^ 10~^ 


to 10~"^ events per galaxy per year (e.g. Magorrian fc 

Tremainel 19991 [Donley et al.|2002[|Wang fc lVlerritt|2L)U4[ 

Kesden|2012[ [Stone fc Metzger|2014[|van Velzen fc Farrar| 

". Although the hrst tidal disruption event (TDE) 


M4| 


candidates were discovered in the X-ray s or UV (e.g 


Komossa et al. 2004 Gezari et al. 2009 Bloom et al. 


recently discovered in the optical (e.g. 

Gezari et al.|2012 

Chornock et al. 

120141 [Holoien et al. 

|2014[ [van VelzenI 

&: h'arrar 2014 

Arcavi et al. 2014|). 

However, the re- 


cent optical IDE 
with theoretical ex 

observations are ditticult to reconcile 

pectations (e.g. Ulmer||1999 IStrubbe 

& Quataert 

to 

o 

o 

Lodato & Rossi 2011), in which the 


The class ical description of a TDE was outlined by 


Rees (1988). In this picture, a star of mass M* and ra- 
dius Hi, approaches a supermassive black hole with mass 
Mbh- The star is disrupted when its parabolic trajec¬ 
tory brings it to a pericenter distance Rp smaller than 
the tidal radius, Rt- The resulting stellar debris has 
a rather narrow distribution of specific angular momen¬ 
tum (all of it close to the specific angular momentum 
of the star), but the distribution of mass with respect 
to specihc binding energy cLM/de is roughly flat from 
e ~ —GMbhR*/Rt to e ~ +GMbhR*/Rt- Because the 
semi-major axis of the most-bound matter is 

a„,i„ ~ (l/2)R^/i?* ^ [MBH/M.y/^RT » Rt (1) 

the orbits of these “tidal streams” are highly eccen¬ 
tric. A uniform distribution of mass per binding en¬ 
ergy implies a mass return rate to the stellar pericenter 
.^fallback cx (t/tp)"^'’'^, where to is th e orbital time of the 
most bound material (|Phinney|1989 ). In this model it is 


further assumed that upon return to pericenter, general 
relativistic apsidal precession causes streams returning 
at different times to shock against each other and dis¬ 
sipate sufficient orbital energy to compress these very 
extended, highly elliptical orbits into approximately cir¬ 
cular orbits with radii ^ 2Rp. The inflow time through 
the accretion disk that then forms is estimated to be 
<C to, so that the accretion rate onto the black hole 
tracks closely the mass return rate of the tidal streams. 
If so, the bolometric light curve should peak at ~ 

2 X s (where A/bh,6.5 = -^bh/10®-^Mq, we 

choose this mass as fiducial because it best matches the 
black hole masses estimated fr om stellar bulge p roper- 
ties for the events observed by Arcavi et al. 2014|) after 
the star is destroyed, reaching a maximum luminosity 
^ 2 X ^ erg/s and then decay oc 


This temporal decay became the hallmark for observa¬ 
tional identification of TDEs, and indeed several TDE 
candidates have been reported as having such a light 
curve. Within this model, the effective temperature of 


the peak would be ~ 4 x K if light from 

the inner rings of the accretion disk reaches distant ob¬ 
servers unimpeded; similarly, its effective radius would 
be ~ lORg ~ 5 X IO^^Mbh. 6.5 cm. Because the outer 
edge of the disk lies at a radius ^ Rt, even the narrow¬ 
est emission lines coming from the accretion flow proper 

would have extremely large line widths, ^ 0 . 2 M 3 'jj g g. 

This simple model faces now serious problems when 
confronted with observations of optical TDE candidates 


(IGezari et aL||2012| 

Chornock et al. 20141 Holoien et aL| 

2014||van Velzen & 

barrar|2014 Arcavi et ai.|2(514|). Up- 


tical light curves of these events show (in rough terms) 


the expected decline, and the rise time agrees 

with the return time expected from a TDE due to a 
^ 5 X 10® Mq SMBH. However, the observed tempera¬ 
ture and bolometric luminosity (2-3 x 10^ K, ~ 10^®- 
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10 "^^ erg/s, taking at face value the published bolometric 
correction) are signihcantly lower than predicted, while 
the inferred black body emission radius is much larger 
(~ 10^^ cm, i.e., ^ 2 x 10^i?gg g). In addition, 
line broadening, which reflects orbital motion, is, at 
^ 10,000 km/s, much smaller than would be expected 
for an accretion disk on the scale of Rt- These results 


tions (as noted also by 

Cuillochon et al.|2014 

Chornock 

et al.||2014 

Arcavi et al. 2014 and others). 


have been criticized by many authors. Some simulations 
find that the heating associated with inter-st ream shocks 


unbin ds a signihcant part of the gas (e.g. Ayal et al. 
20001. If, as is expected for the parameters of many I'UE 


& Ulmer 


1997^ 


events, the peak accretion rate is super-Eddington (Loeb 
'1), the accreti on luminosity itself might 
hi 


power a wind th at expels gas (Ohsuga et ai.|200'5 Dotan 


the observed light curve 

(e-R- 

Strubbe & Quataert 

2009 

Lodato & Rossi 201 1|). 

Lodato et al.j (|2009|) and 

Cuil- 

lochon & Ramirez-Ruiz 

P013) showed that, depending 


essarily constant, leading to further deviations from the 
simple picture. Finally, even if the bolometric light curve 
does follow a decay, the optical light curve should 
not reflect the bolometric on e if the spectrum peak s, as 
it often should, in the EUV ( Lodato fc Rossi|[2011 1 . 


However, there is another aspect of the classicaf TDE 
picture that has thus far escaped much critical review: 
the assumption that the gas circularizes immediately 
upon returning to the vicinity of the black hole, and that 
it does so on the scale of the tidal radius. To circularize 
the returning matter at a radius R < amin, the gas must 
lose an energy ^ GMbh/ 7? per unit mass. Dissipation 
(e.g., in shocks) can help, but only if heating quickly 
leads to true energy loss by radiation that escapes the 
matter. Moreover, in order to dissipate this much energy, 
the shocks must have speeds comparable to the orbital 
speed at radius R. In the earliest work, it was supposed 
that these shocks would take place near the pericenter 
(i.e. at a radius from the black hole ~ Rp and at ap¬ 
proximately the azimuthal angle where the star reaches 
Rp), and w ould be cau sed by strong relativis tic apsidal 
precession (Rees 1988). Reconsideration by Kochanek 


(1994) indicated that any shock near pericenter would be 


much weaker, with the dominant stream convergence due 
to the small angular spread of the tidal streams rather 
than to apsidal precession. In a pseudo-Newtonia n SPH 
simulation coverin g the beginning of mass return, [Ross-| 
wog et al. (2009) also found shocks between the tidal 


streams near apocenter, but because the available kinetic 
energy is least at apocenter, these shocks dissipate rel¬ 
atively weakly. Despite these qualms, it has still been 
generally assumed that somehow the tidal debris would 
quickly “circularize”, i.e., lose enough of its orbital en¬ 
ergy through hydrodynamical processes that it can join 
an accretion disk with a radial extent ^ Rt- 


Recently Shiokawa et al. (2015) took up this question 
using detailed numerical sirnulations. (Because of numer¬ 
ical reasons, these simulations were carried out for a tidal 
disruption of white dwarf by an intermediate mass black 
hole, but the simulations can be scaled to a regular star 


tidally disrupted by a massive black hole.) They found 
that the circularization process is slower than previously 
thought (lasting ~ (5-10)to), produces a flow that is only 
roughly circular ((e) ~ 0.4), and leaves most of the de¬ 
bris at radii nearer a„iin than Rt because the principal 
shocks are located at that scale. By quantitatively defin¬ 
ing “circularization” in terms of progress toward reduc¬ 
ing the eccentricity of the fluid orbits and making the 
flow axisymmetric, they demonstrated that circulariza¬ 
tion may remain incomplete even at the end of the event 
(hence the quotes around the term in the title of this 
paper). Even more importantly, they pointed out that if 
the characteristic radius of the gas is signihcantly greater 
than Rt, inflow could take considerably longer than ~ tpi 
decoupling the time-dependence of the light output from 
the time-dependence of the mass-return rate. 

In this paper, we build upon their analysis of where 
the shocks in this flow are located and the magnitude 
of the heating rate associated with them to point out 
that the very fact they occur at rather larger radius than 
previously expected makes them a strong candi date for 
the origin o f the light seen in optical TDE events (Arcavi 


et al. 
color 


2014). As we will show, the peak luminosity, the 
emperature (and therefore the size of the emitting 


region), and the line widths are all reproduced well by 
this model. 

We use the term “circularizatiorQ’ as a label for our 
outer shocks-powered TDE model in order to distinguish 
it from the commonly assumed accretion-powered model. 
However, we emphasize that in fact the gas orbits remain 
highly elliptical even after passing through these shocks 
because their binding energy remains much smaller than 
that corresponding to a circular orbit with their angular 
momentum. In addition, the gas is partially pressure- 
supported. 

We outline the results of the numerical simulations 
and the apocentric circularization model stemmi ng from 
them in § We summarize the observations (Gezari 
et al. 2012rfChornock et al. 120141 Holoien et al. 


van Velzen fc Farrar||2014[ [Arcavi et al.||2014p ~ an^com- 

pare them to the predictions of this model in § We 
conclude and summarize the implications of these find¬ 
ings in 0 

2 . THE CIRCULARIZATION MODEL 

We begin by briefly summarizing a few basic proper¬ 
ties of tidal disruption events. For a main sequence star, 
the radius R* = RqM\~^, where A4* = M^/Mq and 
^ « 0.2 for 0. 1 < M* < 1, but increases to « 0.4 for 


1 < A4* < 10 (Kippenhahn & Weigert 1994). To account 
for the shape of the star's internal density prohle, we fol¬ 
low Phinney (1989), defining / as the ratio of the gravita¬ 


tional binding energy of the star to GMl /R* and k as the 
apsidal motion constant (determined by the star’s radial 
density profile). The ratio k/f = 0.02 for radiat ive stars, 
but is 0.3 for convective stars (Phinney 1989). These 


ratios are extreme values, as most stars are a mix of ra¬ 
diative and convective regions: in stars less massive than 
the Sun, an outer convective zone surrounds a radiative 

^ Note that [Lodato] (|2012| and [Bogdanovic et al.j ||2014| have 
used the term circularization in the context ot TUtj lightcurves, 
but neither one was referring to the radiation of heat generated in 
shocks at large radius. 
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core; in more massive sta rs, the core is convective whil e 
the envelope is radiative ( Kippenhahn fc Weigert|[l9M ). 
The quantitative balance between the two kinds ot re- 
gions is sensitive to heavy element abundances because 
higher increases opacity per unit mass, enlarging con¬ 
vective zones. This could be a significant effect because 
the stars most likely to be the victims of TDEs are pre¬ 
dominantly drawn from the central regions of galaxies, 
so their heavy element abundan ces may, on average, be 
several times greater than Solar (iSarzi et al.|2005[|R o.ias- 

Arriagada et al.|2014 Gonzalez IJelgado et al. pT^ . We 

will use m the following the geometric mean ot the two 

values, k/f = 0.08, as our fiducial value. 

From the size of the star and a gauge of its internal 
structure (the k/f ratio), we can estimate the tidal radius 


intercepted by matter that has already pass ed through 


Rt 


6.7 X 10^2 

[0.08 ) 


1/6 


M.J ^AfgyggCm. (2) 


In the estimates to follow, we will assume Rp ~ Rt for 
several reasons: the cross section for disruptions is oc 
Rp and furthermore if the stellar loss-cone isn’ t full, the 


shock 1 at least once. Shiokawa et al. (2015) called the 
forward shock, the one in which fresh rnaterial is shocked, 
“shock 3” and the reverse shock, the one acting on gas 
that had already returned, “shock 2”. Over time, shock 2 
divides into two arms, and stream deflection, both by the 
outer shocks and by the increasing pressure in the peri- 
center region, leads to the disappearance of shock 1. The 
mass return rate does roughly follow the classical 
trend, but after 80% of the bound d ebris had returned 
from apocenter, Shiokawa et al.|(2015 ) found that most of 
the mass is placed on only roughly circular orbits at radii 
^ aniin, — (5-8)i?T for their parameters, corresponding 
to amin — 100i?T, when scaled to the tidal disruption 
by a SMBH, which is of interest here. Moreover, a time 
~ (3-10)<o is required for even this degree of “circular¬ 
ization” to be achieved. A fraction of the returning mass 
loses enough angular momentum by shock deflection that 
it is able to accrete within a few toi but the maximum 
a ccretion rate is only ^ O.lx the classical expectation. 


rate of smaller Rp encounters is suppressed (Frank & 


Rees|1976); additionally t he debris energy di s tribution is 


& Ramirez-Ruiz 2013). 


unchanged for Rr, < Rt (Stone et al. 2013 Guillochon 


'i'he semi—major axis of the most-bound material is: 

1/3 


(s) 


Xy^-^M^gggCm. (3) 


The corresponding return time of the most-bound mate¬ 
rial to the pericenter is: 

*0 = 1.8x10= (4) 


and the maximal mass return rate is: 

- 1/2 


Mq ~ —- 

3to 




(5) 


Although Rees (1988) suggested that relativistic ef¬ 
fects cause the apsidal angles for streams with different 
semi-major axes to reach the pericenter with such dif¬ 
ferent directions that their mutual shocks have speeds 
of order the orbits’ maximum orbital speed, quantitative 
study indicates this occurs only if R„/R „ is rather small 
if the initial stellar orbi t is bound ([Bonnerot et al. 


or 


2015 Hayasaki et al. 2015). Similarly, Lense-'l'hirring 


Shiokawa et al. (2015) considered, for numerical rea- 
sons, a 'I'UF ot a white-dwarf by an IMBH. In that case 
amin/Rp ~ 10, and after redistribution of angular mo¬ 
mentum most of the matter stays at a distance com¬ 
parable to Omin. In a MS-SMBH this ratio is of order 
100, and clearly there is not enough angular momen¬ 
tum to keep the matter at a circular orbit there. As 
such, the term “circularization” might be somewhat con¬ 
fusing. However, because energy dissipation is so slow, 
the matter will still remain mostly at a distance com¬ 
parable to Oinin- It will be supported partially by ther¬ 
mal pressure, and it will settle into elliptical orbits. The 
near-apocenter location of the outer shocks still holds 
due to the geometry of elliptical orbits whose apsidal di¬ 
rections are slightly different and because the shocks are 
built upon collision with newly-arriving streams. The re¬ 
orientation of these newly-arriving streams at the loca¬ 
tion of the outer shocks weakens the inner shock and pre¬ 
vents further dissipation, and the streams remain highly 
eccentric in the absence of an additional mechanism for 
energy loss. 

The outer shocks (shocks 2 and 3) dissipate energy 
at a rate comparable to the mass return-rate times the 
returning matter’s orbital kinetic energy. Because the 
density of the matter approaching shock 3 from large 
radius is much greater (by as much as two orders of 
magnitude) than the density of matter that has already 
passed through shock 1 at least once, heatin g shocks 2 


is rather gr eater than in shock 3 (see Fig. 13 Shiokawa 


precession 

also requires relatively small Rn/Rn to be 

significant 

Guillochon & Ramirez-Ruiz 

2015 Hayasaki 


et ah [2015 

. Instead, three separate shock systems (see 

Fig. [IJ), none ot them that strong, combine to transtorm 


et al. (2015)). To order of magnitude accuracy, the outer 
shock heating rate can be estimated by 


the debris motion from highly-elliptical ballistic orbits 
to only moderately ellipt ical orbits significantl y influ- 


E, 


peak 


GMbhMo 




-5/6 


^l/6+5C/2^-l/6 
( 6 ) 


enced by hydrodynamics (Shiokawa et al. 2015). Upon 
passing through the pericenter region, the geometrical 
convergence of strea ms in different orbital planes creates 


the “nozzle shock” (Fvans & Kochanek 
T9M1) dubbed “shock 1” inj^ 
le orb! 


1 989[ Kochanek 
°|l 2015[ ). Intrm- 


Shiokawa et af 
ital apsides between different 


SIC misalignment of the 
stream orbits creates a forward/reverse shock system at 
radii Omin where streams returning for the first time are 


The quality of this estimate can be confirmed by 
scaling the heating rate from the [Shiokawa et al!) 
(2015) simulation, which followed a O.64M0 white dwarf 
(with i?* * = O.OI 24 R 0 ) disrupted by a 5 OOM 0 black 
hole, to conditions appropriate to the fiducial values 
of this paper. We do so assuming that Apeak cx 
M*(GMBH/omin)/Aorb(amin)- The Scaling can then be 
most succinctly accomplished by writing Rt = ^TZR*, 
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with TZ = and 4> = {k /for main se¬ 

quence stars (Phinney 1989); our explicit calculation of 
the white dwarf disruption indicates that $wd — l-l- 
In this notation, the specific energy of the most-bound 
debris e = ~ so that 

Omin = GMBH/(2e) = and Porb(amin) = 


Shock 1 - 




Fig. 1.— Schematic description of the shock system. The inner 
shock 1 takes place near the pericenter, while the outer shocks 2 and 
3 arise nearer the apocenter. Shown are two different times, t ~ 2to 
(top) and t ~ 6to (bottom). Note that at earlier times shock 2 has 
only one branch (top), whereas it splits into two at later times 
(bottom). Sample gas streamlines are shown as black curves with 
arrows; shocks are shown in red. The large black disk is the black 
hole; the thick black curves without arrows are rough indications 
of the boundaries of the flow. Gas-1 denotes returning matter, 
while gas-2 denotes matter that has gone around the black hole at 
least once. The contact discontinuity between gas-1 and gas-2 is 
indicated by a dashed black curve. Separations are not drawn to 
scale in order to emph asize the sequence of events. Adapted from 
[Shiokawa et al.] ( |2015[ |. 



Fig. 2.— Heating rate from scaled to our fiducial main sequence 
parameters as described in text. 


Combining these yields 


E, 


peak 




^BH 


( 7 ) 

The ratio between the predicted heating rate for our fidu¬ 
cial values and that found in the white dwarf simulation 
is then 1.7 x IQ-^ik/ 

assuming ^ = 0.2. Similarly, the timescale ratio is 
8.3 X lO3[(fc//)/O.O8]i/2Mj^/j% 5(M*/M0)O-2. Applying 
these factors to the h eating rate calculated explicitly in 
Shiokawa et al. (2015), we find a heating rate as a func- 
tion of time shown m Figure For our fiducial parame¬ 
ters, it peaks at slightly less man 10^"^ erg s“^, a factor 
of ~ 5 below the order of magnitude estimate of eqn. 
The characteristic time is ~ 21 d, so the peak heating 
rate lasts for several months and then trails off, roughly 
following the late-time decline of the mass-return rate. 
These numbers are more illustrative than general for two 
reasons: some details of TDE stream behavior scale with 
72.; and the simulation assumed a Schwarzschild space- 
time, while an orbital plane inclined relative to a Kerr 
spacetime could introduce further complications (|Guillo-| 


chon & Ramirez-Ruiz||2015 |Hayasaki et al.||2015fT' 

Whether this heat can be radiated quickly depends on 
the local optical depth. Near the apocenter it is 


kM* 


(^) 


-2/3 


( 8 ) 


where we used the Thomson opacity k = 0.34cm^g“^ 
and assumed that only half the star’s mass is bound. 
With such a large scattering optical depth, we assume the 
radiation is well thermalized and escapes with a roughly 
blackbody spectrum. If the vertical scale height 77 ~ i?, 
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the diffusion time tH /c is 


^diff 


rOn 




-1/3 


^ 2 / 3 + 5 ^- 2/3 

This time is, for our fiducial scenario, only a factor of 
three greater than the characteristic timescale of mass- 
return toj but its ratio to to scales with black hole mass 

-7/6 


OC 


M 


RH 


If, as is indicated by the simulation of 


Sh- 


iokawa et al.| (|2015|), H/R ~ 0.5, we expect the lumi- 
nosity to track the heating rate for masses larger than 
^ 3 X 10®Mq, but be delayed with respect to the heat¬ 
ing rate for smaller masses. In the quick diffusion time 
(large black hole mass) limit, the peak radiated bolomet- 
ric luminosity is, up to the efficiency factor, ifpeak- This 
peak luminosity is always smaller than the Eddington 

luminosity LEdd- Although E-peak/Lsdd oc , the 

diffusion time grows slightly faster with decreasing mass. 

For this reason, at lower masses LEdd oc 

so that the maximum L/LEdd is never more than a few 
tenths. Instead, the relatively stronger heating rate when 
the black hole mass is smaller will likely lead to pushing 
matter outward (because the energy liberated is bounded 
by the binding energy, the majority of the mass cannot 
be expelled). 

The pericenter shock also heats the gas, but the opti¬ 
cal depth in this region is so much greater that its con¬ 
tribution to the luminosity is small. Accumulation of 
only 10 “"^ of the infalling matter near Rt results in a 
diffusion time from this region that is longer than the 
accretion time scale (for a = 0 . 1 ). 

The blackbody temperature of the apocenter radiation 
in the rapid diffusion limit is: 


the end of the simulation at f ~ 12 to, its decay is at 
least roughly consistent with which is what would 

be expected if late-time heating were primarily due to 
shocks actingon fres hly- retur ned matter. 

Eqs. Q, M, ( [l^ and ( |TT| provide estimates for 
the values ofnve ^servables (radius, time, Lp^ak, T and 
v) that characterize the emission due to the circulariza¬ 
tion process. To illustrate the likely span of observable 
properties predicted by our model, we computed what we 
would expect from each of three sets of parameters, our 
fiducial set and two others, differing in stellar structure 
(k/f) and black hole mass Mbh- These sample predic¬ 
tions are shown in Table [T] 

We emphasize that the signal discussed here is emitted 
regardless of whatever radiation is produced by other 
processes during the course of the TDE, e.g. the emission 
produced by the accretion process or by outflow-driven 
shocks. For TDEs with Rp/Rg at least a few tens, the 
principal remaining uncertainty is stellar structure; this 
uncertainty is reflected in the factor k/f. On the other 
hand, there may be classes of TDEs to which our analysis 
does not apply. For example, if Rp/Rg < 10, whether 
because Mbh is especially large or Rp is especially small, 
relativistic apsidal precession or, if the black hole rotates, 
Lense-Thirring precession may lead to strong shocks at 
radii nearer (logarithmically) to Rp than to Oniin- In this 
case the circularization radius will be smaller and the 
circularization signal will be brighter and hotter. The 
circularization energy would then be comparable to the 
accretion energy and this signal might blend with the 
accretion signal. 

3. COMPARISON WITH TDE CANDIDATES’ 
OBSERVATIONS 


5.1 X 10' 


Kfht) 
VO.08 J 


-3/8 


( 10 ) 


The typical relative velocity between shocking streams 
at the apocenter region is the Keplerian velocity: 


\ ^min / 


11000 kms 


-1 


0.08 


- 1/6 


Arcavi et al. (|2014 ) describe seven rather similar opti- 
cal IDE candidates. All these candidates (by selection) 
had comparable peak luminosities and light curves with 
similar timescales. Table presents a summa ry of these 
observation s (including in addition P Sl-llaf, Chornock 
et al.|2014 but excluding SDSS J0748, [Wang et al.| 2 l)ll 


M 


If the orbits were circular and the line emission con¬ 
fined to a narrow annulus, the line profile would have a 
pair of peaks separated by 2 i;sin/, for inclination angle i. 
Elliptical motion can alter this separation by a factor of 
order unity and shift its center; a larger radial spread in 
the zone of emission can, to a degree, fill in the minimum 
in the line p rofil e. 

Eqs. Q - (|IT|) agree reasonably wi th the results of the 
simulation by |Shiokawa et al. (2015) when extrapolated 
from their simulated white dwart scenario to our fiducial 
main sequence scenari o. According to the simulation by 
Shiokawa et al. (2015), although shock 1 initially domi- 
nates the dissipation rate, after a time ^ toj shocks 2 and 
3 produce the most heat. The rise of the outer shocks 
marks the beginning of efficient dissipation of orbital en¬ 
ergy at the apocenter region, and it therefore corresponds 
to the rise of the signal that we model. Following the 
peak of the outer shocks’ dissipation, the simulated en¬ 
ergy dissipation rate decreases. Although the heating 
rate falls below the peak rate by only a factor of 2 by 


4 / 2 -i/ 6 T, 4 ^i/ 6 which has no reported Apeak, T and R). Presented in the 
^^^BHtaMe are, for each of these o ptical TDE candid ates, the 
SMBH mass ( as e stimated by Arcavi et al. 2014|using the 


Gadotti 2009 and Haring & Rix|2004 bulge relations), the 
estimated peak bolometric lurninosity, the fitted black- 
body temperature and radius, the width of observed Hell 
or Ha emission lines, and whether the optical light curve 
is consistent with a decay trend. For most of the op¬ 
tical candidates, the SMBH mass is a few times 10®M©. 
For such SMBH masses disrupting a one solar mass main 
sequence star, our model predicts Apeak lO'^^ergs”^, a 
blackbody temperature ^ 4 x lO'^ K, a blackbody radius 
(i.e., amin) ~ 5x10^'^ cm, and aline widths 8000 kms“^. 
These figures agree quite well with the measurements re¬ 
ported in Table 

A unique feature of many optical TDE candidates is 
an observed temperature almost constant in time, in con¬ 
trast with theoretical predictions of accreti on theory. For 


exam ple, the spectral shape of PSl-lOjh (Gezari et al. 
2012) at optical-NUV wavelengths is uncnanged while 
the integrated flux falls by a full order of magnitude, and 


the ASASSN-14ae (Holoien et al. 2014) fitted tempera¬ 
ture remains ~ 2 x 10'‘ K during a fall of one and a half 









































6 


TABLE 1 

Model predictions for TDEs with a rise time of a month 


Stellar 

structure 

k/f 

A^bh 

(IO 6 M 0 ) 

Apeak 

( 10 ^^ergs“^) 

Tbb 
(10^ K) 

A?bb 

( 10 ^® cm) 

Line width 

(kms“^) 

radiative 

0.02 

10 

130 

5.6 

0.44 

17000 

convective 

0.3 

1 

20 

4.8 

0.23 

7500 

fiducial 

0.08 

3 

50 

5.1 

0.31 

11000 


The first prototype matches predominantly radiative main sequence stars, those with 
A4* ~ 1. The second prototype is for mostly convective stars, which tend to be 
stars of either very low or very high mass. Our fiducial case, presented in the third 
line, is defined by the geometrical means of the extreme values for the k/f factor 
and the black hole mass. Apeak presented here and the corresponding temperature 
(see Eqs. [^and |10[ l assume 100% efficiency of the shocks in converting gravitational 
energy to observea lu minosity. The actual efficiency is lower, e.g. a comparison to 
[Shiokawa et al.|2015| suggests luminosities lower by a factor of p:: 5 and temperatures 
lower by a factor ot ^1.5 (see text). 


TABLE 2 

Observed properties of optical TDE candidates 


Event 

A4bh 

Apeak 

Pbb 

A?bb 

Line width 


(10® M 0 ) 

(10^®ergs“'^) 

(10^ K) 

(10^® cm) 

(kms“^) 

PSl-lOjfil 


^ 22 

^ 3 

0.6 

9000 ± 700 

PSl-llaf2 

8± 2 

8.5 ± 0.2 

1.90 ± 0.07 

~ 1.2 


PTFOOge® 

0.00_Q gg 

Q t: + 50 

ot)_40 

3.1 ± 0.3 

1.14 ± 0.2 

10070 ± 670 



^+5.3 (a) 

3.3 

2.2 ±0.3 

0.59l°;l® 


SDSS TDE2'‘ 

Qcr c-l-55.3 

> 4.l(‘’) 

1 Q9-I-O.O7 

t-o^_0.06 


-80001=1 

ASASSN14ae® 

2 45"^^’^^ 
^•^^-0.74 

8.2 ± 0.5 

2.2 ± 0.1 

0.7 ± 0.03 

17000-80001=) 

PTFOOaxc® 

^•oy_o.64 

Q+3.3(d) 

^•^-1.4 

-^•-^^-0.17 


11890 ± 220 

PTFOOdjl® 

0 ir7+9.97 

-2.96 

^ „+23.1(e) 

' -10.4 

9 7-I-O.7 
'-0.5 

0.58l°;24 

6530 ± 350 


fl) [Gezari et alT 


. _ . (2f iChornock et al.|2014| (3) lArcavi et al.|2014| (4) Ivan Velzen &:! 

|Farr'ar|^Ul4 T57l tloloien et ai.|2Ui4| “ 

Lor ttie nrst four events tne autnors reported a reasonable agreement between the ob¬ 
served fight curve and a ~ decline. ASASSN-14ae does not show such a decline 

while the last two are not clear. 

(a) PTF09ge: First fine, values we fitted from the peak band luminosities kindly supplied 
bv lair Arcavi: second line, oublished values for L. Trr and i?RR at 19 davs before peak 
dArcavi et al.|20l4. 

\bj I'lie l^peak values presented are bolometric, except for SDSS TDE2, where it is only 
the g band. The bolometric peak for this event is not published, but is probably similar 
to or higher than PSl-lOjh. 

fc) For SDSS TDE2 and ASASSN-Mae, the values are from the discovery papers (|vaR| 

IVelzen &; Farrar|2014||Holoien et al.|2014|l. These should be comp ared with 3440 di llUU 
and ^bUU ± LYb (respectively), reported by [Arcavi et al.| | |2014| t. For ASASSN-14ae, 
17000 8000 means 17000kms~^ at the luminosity peak and a decrease to 8000kms“^ 

three months later. 

fd) A valu e estimated using the published values of T and H (determined by [Arcavij 
jet al.|2014| using a fit to the continuum at around the peak of the fight curve). Note that 
tnis is coniparable to the r-band luminosity and thus should be taken as a lower limit. 

fe) A valu e estimated using the published values of T and H (determined by [Arcavij 
[et al.|2014'| using a fit to the continuum at around the peak of the light curve). 


orders of niagnitude in integrated flux. When the tem¬ 
perature is constant, the bolometric correction should be 
likewise constant. A decay in the luminosity combined 
with a constant T then implies that the radius is decreas¬ 
ing with time. A possible explanation is a gradual shift 
inward of shock 2 as debris that already passed through 
pericenter gradually settles into less eccentric orbits. 

One can easily see the general agreement between the 
observations reported in Table 3 and t he m o del predic¬ 
tions sjrmmarized in Eqs. (|3l, (4|, ([^, ([To|, ([TTl) and in 
Table The values of the estimate luminosities (and 


corresponding temperatures) are slightly larger than the 
observed values. This is reasonable given that after all 
they are only and order of magnitude estimates which 
assume 100% efficiency. In deed a comparison wit h the 
scaled luminosities based on Shiokawa et al. (2015) show 
indeed luminosities that are smaller by about a factor of 
5. Although the observed values do not exactly match 
the values expected for our fiducial scenario, we did not 
attempt to fit the model parameters to the observations. 
Both the model and the data (N.B. the absence of extinc¬ 
tion corrections and the discrepancies between different 
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observations noted in Table are too crude for that. 
The rise time to peak, roughly reflecting to, is about a 
month for all events with a well observed rise phase. This 
may reflect a selection effect, as many of the events were 
found in data of surveys with cadences optimized for de¬ 
tection of supernovae. 

For any particular stellar structure, i.e., choice of k/f, 
any pair of the four observables Tpeak, T, v and to could 
in principle be used to solve for Mbh and M*, thus over¬ 
constraining these values. However, varying M* within 
a factor of a few, reasonable for main sequence stars, has 
only a marginal effect on the observables, which depend 
very weakly on the mass of the tidally disrupted star. 
Therefore, given a value of fc//, the observables’ values 
are mostly determined by the SMBH mass. For example, 
for our fiducial value of fc//, a rise time of to 1 month 
implies Mbh,6 .5 ~ 1- In addition, both larger k/f and 
larger Mbh diminish the luminosity and temperature, 
but increase the line-width. Consequently, to match a 
specific range of observables, a smaller fc// (a star more 
of whose volume is radiative) demands a larger Mbh, 
and vice versa. However, the degree of sensitivity to 
these parameters varies; in particular, the line width and 
peak luminosity depend more strongly on Mbh than the 
temperature. Thus, there are distinctions between the 
three cases listed in Table [TJ 


4. DISCUSSION 

We suggest here that the observed optical light from 
TDEs arises from emission on radial scales ~ a„iin, and 
that the energy source is the orbital energy dissipated 
by shocks at that distance from the central SMBH. Al¬ 
though the energy dissipated in this process is a small 
fraction (~ 10“^) of the total energy available from ac¬ 
cretion onto the black hole, it is actually of the right 
scale to power the observed signal. The observed tem¬ 
peratures, emission radii and line widths are all in agree¬ 
ment with the expected characteristics of this emission. 
The observed rate of decline is also crudely consistent 
with the decline in heating rate due to these shocks. All 
this, plus the fact that radiation from the region of tidal 
stream shocks must occur in any event, are significant 
advantages. 

At the same time, however, these successes point to 
a puzzle. Our model says nothing about what happens 
to the ~ 10®^ erg one might expect to be released when 
the debris is accreted onto the black hole. Where does 
this energy go? In fact, this question is raised equally 
strongly in regard to models in which all the radiation 
stems directly from accretion. Posed in that context, it 
becomes: why do we see only 1% of the expected energy? 
The only difference between the nature of this question 
as it applies to our model and as it applies to the more 
conventional models is that in our case the suppression 
factor could be anything smaller than ^ O.OI; in the 
conventional models, it must be consistently ~ 0.01. 

A number of speculative answers to this question ex¬ 
ist. One possibility is that the mass return rate can of¬ 
ten be super-Eddington, suggesting strong photon trap¬ 
ping in the deepest part of the gravitati onal well, where 
most of the accretion ene rgy is released (Begelman|1979 
Abramowicz et al.|1988). This o ption raises seve r al dit- 
hculties. On the theoretical side, Shiokawa et al. (20151 
found that the peak accretion rate onto the black hole 


is reduced by an order of magnitude relative to conven¬ 
tional models, restricting the parameter space of super- 
Eddington accretion to lower-mass black holes. In ad¬ 
dition, recent simulational work on super-Eddington ac¬ 
cretion has shown that magnetic buoyancy effects may 
provide an end-run around photon-trapping that p ermits 


the emergin g luminosity to match the heating rate (Jiang 
et al.|201^ . On the observational side, if strong photon- 
trapping do 


trapping does occur, it would create a light cur ve with 


a lengthy peri od of nearly constant luminosity (Krolik 


& Piran 20121, and none of the events discussed here 


shows such behavior. Indeed, the luminosity begins to 
decline within a month or two of the peak. If photon¬ 
trapping explained the small radiated energy, this com¬ 
mon behavior would require fine-tuning, for all observed 
events would have to be just slightly super-Eddington. 
Note, however, that the peak mass-return ra te in Ed- 


dington units is only ~ ^ (Shiokawa 


et al. 2015); there may be a population of events with 
somewhat larger masses for which the mass-return rate 
is super-Eddington for only a short period of time. An¬ 
other option is that when the mass accretion rate is 
super-Eddington, 99% or more of the accretion power 
is put into ki netic energy of a low-mass outflow ra ther 
than photons ( Ayal et Sl]| 2000 Ohsuga et al. |2005D . In 
the super-E ddington accretion simulations of p^dowsk” 


et al. (2014), the effective radiative efficiency in rest-mass 


units at an accretion rate lOOMEdd (defined with respect 
to an assumed radiative efficiency of 0.057) falling into a 
non-spinning black hole was indeed depressed by about a 
factor of 100. However, when they examined such a flow 
onto a black hole with spin parameter 0.9, they found 
a suppression factor of only about 20 and a substantial 
jet power. Their simulation also raises another problem 
with stronger implications for this option: they predict a 
very strong anisotropy of the emitted photon luminosity. 
One would then expect that surveys would strongly favor 
discovering systems in which the apparent radiative effi¬ 
ciency is much greater, and their spectra would be much 
hotter than those actually seen. 

A third is that 99% or more of the accretion flow is 
blown away before it ever comes near the black hole. 
In one version of this suggestion, the wind prevents 
the accretion rate ont o the black hole from si gnificantly 
exceeding Eddington (Poutanen et al. 2007). In that 
case, it suffers from the same lightcurve problem as 
photon-trapping, for both amount to mechanisms that 
regulate the accretion onto the black hole to a steady 
near-Eddington level, for any value of the accretion rate 
greater than Eddington. It has also been suggested that 
the mass-loss might be regulated by atomic line opac¬ 


ity (Laor & Davis 2014 Miller 20I5|). This version has 
the desirable property or predicting tne luminosity to de¬ 
crease as the accretion rate diminishes, but not as rapidly 
as the magnitude of the accretion rate in the outer disk 
falls. On the other hand, although there are significant 
intrinsic uncertainties in this model’s prediction of the 
luminosity, in rough terms it predicts a total radiated 
energy only a factor ~ 10 smaller than the conventional 
prediction because nearly all the accreting matter reaches 
radii ~ 100i?g 

A fourth imaginable answer is that the additional en¬ 
ergy is radiated away in the EUV/FUV, or non-thermally 














































in a different channel that has not yet been observed. 
Such a situation would be expected if, when accretion 
onto the black hole finally begins, the light generated in 
the inner disk avoids reprocessing on its way out. An¬ 
other version of this explanation is that the actual ex¬ 
tinction corrections are much larger than currently sup¬ 
posed. Without further definition, it is hard to evaluate 
this class of solution. 

Finally, we comment that the lengthened timescale for 
any particular com bination of stella r mass and black hole 
mass predicted by Shiokawa et al. (20151 bears on this 
question only in the sense that it dirnimshes the range of 
parameters for which super-Eddington accretion would 
be expected. The total accretion energy predicted would 
be the same as any other model in which roughly half of 
the disrupted star is ultimately accreted onto the black 
hole. 


We close with a few thoughts about observations that 
could provide tests for our model. Our predicted peak 
luminosity depends very weakly on the black hole’s mass, 
but the temperature decreases and the expected line- 
width increases when Mbh increases. Thus, a sample of 
events with well-determined black hole masses, peak lu¬ 
minosities, emission line widths, and characteristic tem¬ 
peratures would provide a strong statistical test. 
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